Nonlinear AC resistivity in s-wave and d-wave disordered granular superconductors 
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We model s-wave and d-wave disordered granular superconductors with a three-dimensional lattice 
of randomly distributed Josephson junctions with finite self-inductance. The nonlinear ac resistiv- 
ity p2 of these systems was calculated using Langevin dynamical equations. The current amplitude 
dependence of pi at the peak position is found to be a power law characterized by exponent a. The 
later is not universal but depends on the self-inductance and current regimes. In the weak current 
regime a is independent of the self-inductance and a = 0.5±0.1 for both of s- and d-wave materials. 
In the strong current regime the values of a depend on the screening. We find a ~ 1 for some inter- 
val of inductance which agrees with the experimental finding for d-wave ceramic superconductors. 
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The symmetry of the superconducting pairing func- 
tion has been of great interest lately. The gap of con- 
ventional superconductors has s-wave symmetry whereas 
there is now good evidence that the superconducting gap 
of the high-T c cuprates has d-wave symmetry pj] . Gran- 
ular superconductors are usually described as a random 
network of superconducting grains coupled by Josephson 
links Pll- In high-T c ceramics, depending on the rela- 
tive orientation of the d-wave superconducting grains, it 
is possible to have weak links with negative Josephson 
coupling Q, which are called 7r-junctions. The existence 
of these 7r- junctions may cause, e.g., the paramagnetic 
Meissner effect 0] observed at low magnetic fields ||. 

Recently, Kawamura || proposed that a novel thermo- 
dynamic phase may occur in zero external magnetic field 
in unconventional superconductors. This phase is charac- 
terized by a broken time-reversal symmetry and is called 
chiral glass phase. The frustration effect due to the ran- 
dom distribution of tt junctions leads to a glass state of 
quenched-in "chiralities" , which are local loop supercur- 
rents circulating over grains and carrying a half-quantum 
of flux. 0| Evidence for the transition to chiral glass has 
been seen from experimental studies of the nonlinear ac 
magnetic susceptibility f| , the dynamic scaling || and 
the aging phenomenon fllQfl . The susceptibility measure- 
ments of Ishida et al. [|ll| do not, however, support the 
existence of the chiral glass. 

In order to further probe existence of the chiral glass 
phase Yamao et al. |12| have measured the ac linear re- 
sistivity po and the nonlinear resistivity p 2 of ceramic 
superconductor YL^CuiOs- P2 is defined as the third 
coefficient of the expansion of the voltage V(t) in terms 
of the external current I ext (t): 



V = p I e 
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P2 = -4V^// 3 
V(t) sm(3ut)d{ut) 



(2) 



When the sample is driven by an ac current I ex t(t) 
Iosm(ujt), one can relate p 2 to third harmonics 
the following way 



Yamao et al. have made two key observations. First, 
since the linear resistivity does not vanish at the peak 
position of p 2 they identify the transition as a transition 
to the chiral glass phase. Second interesting observation 
is the power law dependence of \V^(T p )/I ) 3 \ (or of p 2 ) 
at its maximum position T p on Iq; l^'a>(^p)/^o I ~ AT"- 
The experimental value of the power law exponent was 
a w 1.1. Using the XY-like model for <i-wave supercon- 
ductors Li and Dominguez p3| were able to reproduce 
the experimental results of Yamao et al. Ji2| qualitatively. 
The quantitative agreement was, however, poor and the 
role of inductance was not explored. Namely, a was com- 
puted only for one value of dimcnsionlcss inductance L=l 
and with large error bars [13]: a = 1.1 ± 0.6. 

The goal of this paper is twofold. First, we calculate 
a with high accuracy for both of s- and d-wave systems 
using the Langevin equations for the XY-like model with 
screening. Second, we try to answer the question if it is 
possible to discriminate between s- and d- pairing sym- 
metry by measurements of a. We show that there are 
two distinct current regimes for a. In the weak current 
regime (WCR) (small Io) this exponent does not depend 
on the inductance and a = 0.50±0.1 for s- and d-wave ce- 
ramics. In the strong current regime (SCR) a depends on 
the screening. For small L we obtain ad-wave > ot s - wave , 
possibly because in the weak screening limit the energy 
landscape of the d-wave case is more rugged than the s- 
wave case. As the self-inductance increases the number 
of energy local minima gets smaller JLlj and the behavior 
of the two systems becomes more similar, with the values 
of a being almost the same. For the d-wave system in the 
SCR and with 1 < L < 5 we find a rj 1.0 which agrees 
with the experimental value p"2[ . 

We consider the following "coarse grained" Hamilto- 
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where 0, is the phase of the condensate of the grain at 
the i-th site of a simple cubic lattice, Jij denotes the 
Josephson coupling between the i-th and j-th grains, L 
is the self-inductance of a loop (an elementary plaque- 
tte) , while the mutual inductance between different loops 
is neglected. The first sum is taken over all nearest- 
neighbor pairs and the second sum is taken over all el- 
ementary plaquettes on the lattice. Fluctuating vari- 
ables to be summed over are the phase variables, 0$, at 
each site and the gauge variables, A^j = ^ J? A(r)dr, 

at each link. $ p = ^ S<y > Aij ^ s the total magnetic 
flux threading through the p-th plaquette, and <f>o de- 
notes the flux quantum. The effect of screening currents 
inside grains is not considered explicitly, since for large 
length scales they simply lead to a hamiltonian TL with 
an effective self- inductance L |l7| . 

For the <i-wave superconductors we assume Jij to be 
an independent random variable taking the values J or 
—J with equal probability (± J or bimodal distribution), 
each representing and ir junctions. For the s-wave su- 
percondutors Jij is always positive but distributed uni- 
formly between and 2 J. It should be noted that model 
(3) with uniform couplings was first studied by Dasgupta 
and Halperin pl|. Random 7r-junction models (in which 
Jij is allowed to take negative values with certain prob- 
ability) have also been adequate to explain several phe- 
nomena observed in high-T c superconductors such as the 
anomalous microwave absorption, |T^,|l4| the compensa- 
tion effect |^0j , the effect of applied electric fields in the 
apparent critical current (2lJ and the aging effect. 22 



In order to study transport properties, we use the re- 
sistively shunted junction model. j| Then in addition to 
the Josephson current one has the contribution of a dis- 
sipative ohmic current due to an intergrain resistance R 
and the Langevin noise current. We have redefined nota- 
tion: the site of each grain is at position n = (n x , n yi n z ) 
(i.e. i = n); the lattice directions are /j, = x, y,z; the 
link variables are between sites n and n + fi (i.e. link 
ij = link n, fi) ; and the plaquettes p are defined by the 
site n and the normal direction fi (i.e plaquette p = pla- 
quette n, /x, for example the plaquette n, z is centered at 
position n + (x + y)/2). Then the gauge invariant phase 
differences OJn) = A+0(n) — A M (n) obey the following 
equations fiEa]: 

h d0 M (n) 2e . 

^R-^tT = - ¥ W^(n)-<W^ t 

-A A ; [A+0„(n) - A+0„(n)] -^(n,*) , 
MrM^Ki')} = _V<W^-0 , (4) 



where rj^(n,t) is the Langevin noise current. The for- 
ward difference operator is Atfl v (n) = 9 V (n + fx) — 0„(n) 
and the backward operator A^#„(n) = ^(n) — u (n — y). 
In what follows we will consider currents normalized by 
I j = 2eJ/h, time by r = (j) /2TTjR, voltages by RIj, 
temperature by J/ks and inductance by <fio/2irJ. Free 
boundary conditions and numerical integration are im- 
plemented in the same way as in [|f3 15 1. Depending on 
values of Iq and lo the number of samples used for the 
disorder- averaging ranges between 5 and 800. The num- 
ber of integration steps is chosen to be 10 5 — 5 x 10 5 . 
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FIG. 1. a) Upper panel: the temperature dependence of 
Vsu/Iq for the s-wave system. System size I = 8,L = 1 and 
Jo = 0.1. The open triangles, squares and hexagons corre- 
spond to u = 0.001,0.0005 and 0.0002. The peak is located 
at T p — 1.4. b) lower panel: the same as in upper panel but 
uj = 0.001. The open triangles, squares and hexagons corre- 
spond to Io — 0.2, 0.1 and 0.05. The results are averaged over 
15 - 40 samples. 

The temperature dependence of the nonlinear resistiv- 
ity pi of the s-wave system for Io = 0.1 and for different 
values of u> is shown in upper panel of Fig. 1. Similar 
to the d-wave case [|l3| , there is no visible dependence on 
lo. As seen in lower pannel, as Iq decreases peak values 
of p 2 tends to diverge. For L = 1 the peak is located at 
T p = 1.4 and it coincides with the metal - superconduc- 
tor transition at which thermodynamic quantities diverge 
and the linear resistivity po vanishes. It should be noted 
that our disordered s-wave model is different from the 
gauge glass model p| (in the later case the screening 
spoils the transition to the superconducting state). Fig. 
2 shows the Io dependence of max IV^/TqI of the s-wave 
samples (L=l). Clearly, we have two distinct regimes 
for small and large currents. In the WCR (Iq < 0.1) 
a = 0.50 ± 0.04 and a = 0.51 ± 0.03 for I = 8 and 
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I = 12, respectively. In the second regime we obtain 
a = 1.0 ± 0.05 and a = 1.07 ± 0.02 for I = 8 and I = 12, 
respectively. Since within the error bars the finite system 
size effect is negligible, we will consider only the system 
size / = 8. 




In I 

FIG. 2. The current dependence of max\V3 U / Iq \ for s-wave 
superconductors. We choose uj = 0.001 and L = 1. In the 
WCR a = 0.50 ± 0.04 and 0.51 ± 0.03 for the system size 
/ = 8 and 12, respectively. In the SCR a = 1.0 ± 0.05 and 
a = 1.07 ±0.02 for I = 8 and I = 12, respectively. The results 
are averaged over 5 - 800 samples. 

Fig. 3 shows the dependence of max |V^,/£g| on I 
for the d-wave case (1 = 8 and u> = 0.001). In the 
weak current part one has a = 0.51 ± 0.03, 0.45 ± 0.05, 
0.48 ± 0.05 and 0.43 ± 0.06 for £=0.1, 1, 10 and 20, re- 
spectively. Clearly, within error bars a is not sensitive 
to the screening. In the SCR it becomes dependent on 
L: a = 1.8 ± 0.16, 1.56 ± 0.17, 0.97 ± 0.02 and 0.60 ± 0.02 
for £,=0.1, 1, 10 and 20, respectively. Fig. 4 shows the 
results obtained in the SCR for s- and d-wave systems 
with different values of £. The power law region of the 
rf-wave case is sensitive to the screening and is narrower 
than its s-wave counterpart. 

The dependence of a on £ in the SCR is shown in Fig. 
5. Such a dependence may be understood taking into 
account the interplay between the thermal fluctuations 
and the rugged energy landscape. In the weak screening 
limit the later plays an important role and a of the d- 
wave system is bigger than that for the s-wave one. As 
£ increases the thermal fluctuations take over and the 
opposite situation would happen. The pronounced dif- 
ference between two types of symmetry is seen only in 
the weak screening region. 

It is tempting to interpret the two regimes for a as the 
WCR corresponding to the critical regime for P2(T c ,Iq) 
(since £o — » 0) and the SCR corresponding to a mean- 



field regime (away from criticality) . If there is a contin- 
uous phase transition at a critical temperature T c = T p , 

then current-voltage scaling Q predicts that V ~ £ 3 = T 
at T c , with z the dynamical exponent. Therefore, the 

4±i-3 

non-linear resistivity should be P2(T C ) ~ £ d_1 , and 
thus the expected WCR value is a — (5 — z)/2 in d = 3. 
This predicts that a peak in P2{T) at T c is possible if 
z < 5 (i.e. a > 0). In our case, we obtain a w 0.5 and 
therefore z 4 for the disordered s-wave transition. 
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FIG. 3. The current dependence of max IV^/Iq \ for d-wave 
system. We choose the system size I — 8, u — 0.001 and L= 
0.1, 1, 10 and 20 (its values are shown next to the curves). 
For each inductance one has two distinct current regimes. The 
results are averaged over 10 - 800 samples. 

In the experiment of Ref. jl2| the temperature T p 
is merely an intergrain ordering transition temperature 
above which the thermoremanent magnetization disap- 
pears. In the previous simulations of Jl3| for the <i-wave 
system, T p is the temperature where there is an onset 
of positive magnetization, i.e. the paramagnetic Meiss- 
ner effect starts to be observed, but it does not seem to 
correspond to a phase transition. The chiral glass phase 
transition temperature T cg is found at a lower tempera- 



ture, Tea < T p (for £ = 1, e.g., T C; 



0.29 



Kawa- 



mura p5[ | has found that z w 6 > 5 for the chiral glass 
transition, and thus no peak in P2(T) is expected for this 
transition according to the scaling argument. Therefore, 
the peak measured by Yamao et. al. may not correspond 
to the chiral glass transition, but to the crossover we find 
at T p for the d-wave case. 

In order to compare our results with experiments we 
first show that Yamao et. al. |12| performed measure- 
ments in the SCR. Since real current is / 
J ~ 10 2 K and I - 10" 1 we have I - 10" 
the other hand, the current used in experiments I ~ 10 
mA suggests that the experiments were performed in the 



_ 2eJ t 

2 mA. On 



3 



SCR. As seen from Fig. 5, the value of a in the SCR 
for 1 < L < 5 coincides with the experimental value 
p2| . This interval of inductance is realistic for ceram- 
ics [^6| because typical values of L are bigger than 3. An 
accurate comparison between theory and experiments re- 
quires, however, the knowledge of L which is not known 
for the compound of YBa2Cu40s studied in Ref. E^]. 

-1 



x 



-2 

-3 
-1 



SCR 



0.5 



15 



d— wave 



s— wave 



..... 0.5 



15 



_L 



-1.5 



-1 



-0.5 



Ln L 



FIG. 4. The current and self-inductance dependence of 
max\V3 U /Io\ for d- and s-wave systems in the SCR for L= 
0.5, 5 , and 15 (they are shown next to the curves). We 
choose the system size I = 8 and ui = 0.001. The results are 
averaged over 5-10 samples. 
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FIG. 5. Dependence of a on L in the SCR for s- and d-wave 
systems. 

In conclusion, we have calculated the non-linear ac re- 



sistivity exponent a for s and d-wave granular super- 
conductors with high accuracy. Our results reveal two 
distinct current regimes. In the WCR a is independent 
of the screening strength and of types of pairing symme- 
try. In the opposite case this exponent depends on L. 
A difference between s- and d-wave symmetries in the 
nonlinear resistivity can only be found in samples with 
weak screening. The agreement between simulation and 
experimental results is possible for some interval of L. 
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